Dephasing of electrons in the Aharonov-Bohm interferometer with a single-molecular 
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Abstract 

Phase relaxation of electrons transferring through an electromechanical transistor is studied using 
the Aharonov-Bohm interferometer. With the approach of quantum master equation, the phase 
properties of an electron are numerically analyzed based on the interference fringes. Coherence 
of electron is partially destroyed by its scattering on excited levels of the local nanomechanical 
oscillator. Transmission amplitudes with respect to two adjacent mechanical vibrational levels have 
a phase difference of 7r. The character of phase shift by n depends on the oscillator frequency 
only and is robust for the wide range variance of the applied voltage, the tunneling length and the 
damping rate of the mechanical oscillator. 

PACS numbers: 85.35.Ds, 85.85+j, 63.20.kd, 73.23.Hk 
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1. Introduction 

Discrete quantum states and coherent electronic trans- 
port are two important properties of mesoscopic conduc- 
tors. In a recent experiment of a single- molecule Cqq 
transistor, current steps because of the quantum behav- 
ior of nanomechanical motion was observed TJ. In the se- 
quent years, the electron transport through the mechani- 
cal vibrational junctions becomes a subject of much top- 
ical interests. Such systems were fabricated with a gold 
nanoparticle oscillator quite recently More fascinat- 
ing features of the electromechanical systems have been 
predicted with theoretical approaches, such as negative 
differential conductance 0, 0] , shuttling effect 0, @| , su- 
per and sub-Poissonian Fano factor [7H10( , and spintronic 
transport pTl4l3j . etc. 

By applying the Aharonov-Bohm (AB) interferometer 
with a quantum dot(QD) embedded in one arm, coher- 
ence character of electron tunneling through the QD is 
studied in experiments [l4l - fl6l |. These experiments show 
that a fixed QD supports coherent transport and causes 
a phase-shift to an electron. When a QD is allowed 
to mechanically oscillate around its equilibrium point, 
an electron transferring through the dot would be ac- 
companied by absorbtion or emission of phonons ran- 
domly. Phase property in the mechanical vibration as- 
sisted electron tunneling is still an open interesting ques- 
tion. The vibrational motion in the system can be mod- 
eled as a monochromatic oscillator with a good approx- 
imation. It is different from thermal fluctuating bosonic 
baths which cause decoherence to local electronic states 
of charge [13, EH and spin ■ As recently reported, 

a single vibrational mode of QD-array enhances elec- 
tron transport and partially preserves its phase informa- 
tion (2l| . It is worth mentioning that coherent transport 
of electrons in QDs is also sensitive to spin flip, electron- 
electron interaction and external detectors [22h30| . 

In this work, we shall investigate dephasing induced by 
the electromechanical vibration in the single-molecular 
transistor. It is implemented by embedding a harmoni- 



cally movable QD in one (target) arm of the AB interfer- 
ometer and locating a fixed QD in the other (reference) 
arm. The reason of using the QD in the reference arm 
is that phase shift corresponding to each discrete level 
of the mechanical oscillator can be observed by chang- 
ing the gate voltage of the reference QD. The previous 
research close to our issue of interest is the which-path 
detector of charge by using a cantilever [3l|, [Hj]. This 
detector is based on dot-cantilever coupling which causes 
remarkable dephasing to electrons. In their model, the 
dot-lead coupling does not depend on oscillator position, 
whereas the position dependance is considered in our sys- 
tem since it is significant for the electromechanical shut- 
tle junction 0, Q. In our previous report, we derived a 
full quantum mechanical master equation to describe the 
electromechanical system [l(|. In the equation we con- 
sidered both diagonal and off-diagonal couplings between 
electron tunneling and vibrational mode, and revealed 
that the off-diagonal couplings have important contribu- 
tion to the electronic current. Here, we shall develop this 
approach to describe the influence of the electromechan- 
ical system to the AB interference. To obtain knowledge 
of the phase properties of electron scattering on the vibra- 
tional junction, interference transmission as a function of 
an external magnetic flux through the AB ring will be an- 
alyzed with respect to various parameters. This article 
will discuss in the following that coherence of electron 
in the AB ring is suppressed due to excitation of the 
vibrational mode in the transport process. The suppres- 
sion becomes more serious when one increases the bias 
voltage or decreases the oscillator damping rate. It is 
shown that by moving the energy level of the reference 
QD, global phase shift in the transmission amplitude can 
be observed from the change of interference fringe. The 
transmission amplitudes corresponding to two neighbor- 
ing resonant levels of the molecular junction are off phase 
by 7r. This phase difference causes destructive interfer- 
ence to propagating waves and destroy the coherence of 
electron. 

This paper is organized as follows: In section 2, we de- 
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rive the master equation for the model of the AB inter- 
ferometer with a single-molecular transistor in one arm. 
In section 3, the method of our numerical calculation is 
introduced. In section 4, we present solutions for the var- 
ious parameters. Influence of the electromechanical sys- 
tem to the interference fringe is analyzed. In addition, 
the phase shifts of electron propagating waves through 
the resonating QD will be illustrated. Finally, conclu- 
sions are given in section 5. 

2. Model and equation of motion 

The schematic structure of our AB interferometer is 
illustrated in figure [T] It contains two single-level QDs 
coupled to two electronic leads in parallel. One of them 
with energy s\ (QD1) is fixed in the upper arm and the 
other with energy £2 (QD2) is localized in the lower arm. 
The two arms and the electrodes enclose a magnetic flux 
$ passing through the loop-plane. The QD2 is assumed 
to be bounded in a harmonic potential, which consists 
the electromechanical shuttle junction. The QD1 in the 
upper arm provides reference path. We consider both 
inter and intra-dot Coulomb blockade limits in order to 
make sure that electrons propagate through the two-path 
interferometer one by one. Spin degree of freedom is 
not involved in our approach. The Hamiltonian can be 
written in the form of 

H = Hi ea d s + Hdots + Hmech + H tun , (1) 

where 

Hl ea ds = ^ £ykdy k d y k (2) 

k;y—l,r 

describes the noninteracting electrons in the left (y — 
I) and right (y=r) leads. d yk and d y k are creation and 
annihilation operators of electrons with momentum k and 
energy In Hamiltonian 

Hdots = eicjci + (e 2 - + a))c\c 2 , (3) 

c\(ci) is the creation (annihilation) operator of QDi 
(i=l,2). Here, the last term is the contribution from 
the work of an electric field on charged QD2. The cou- 
pling coefficient is given by = eVxo/hd with the bias 
voltage V and effective distance d between the two elec- 
trodes, e is the absolute value of the electron charge a nd 
xq is the zero point position uncertainty \J h/2mujQ of 
the oscillator with frequency ujq and effective mass m. 
The nanomechanical vibration is treated in the quantum 
regime as 

Hmech = hw a)a + ftu k blb k + ^ h(gb\a + h.c). (4) 
k k 

a (a^) and bk (b k ) are annihilation (creation) operators 
for the vibrational mode and its thermal bath, respec- 
tively. u>k denotes frequency of mode k in thermal bath 
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FIG. 1: Two single-level quantum dots connect to two leads 
in parallel, which enclose a magnetic flux $ for Aharonov- 
Bohm interference. The upper dot is fixed and the lower 
dot is bounded by an harmonic potential and movable in the 
horizontal direction. 

which is coupled to the oscillator with a coefficient g. 
Tunneling through the two QDs is represented by 

Htun (T 2y e- lS ^/V^ at+a >4 fe c 2 + h.c.) 

k;y=l.r 

+h J2 (T lv e iS »*' 4 dl kCl + h.c.) (5) 

k;y=l,r 

The tunneling coefficients between the two leads and 
QD1 are given by T ly e iS ^l A {S l{r) = -1 (+1)) and 
its complex conjugate, where the phase <f> is related to 
the magnetic flux <j) = 27r$/$o with the flux quantum 
$0 = h/e. The coupling coefficient with respect to QD2 
is written as T2 V e~ lSy ' t '/ i e Sya ^ a +a ) which exponentially 
depends on the position of the oscillator. The parameter 
a is inverse ratio of the tunneling length A as a = xq/X. 

We describe the interference process by the quantum 
master equation. It is extended from an equation in our 
earlier publication where a more detail derivation can be 
found [10| • The state of total configuration is written by 
the density matrix pt (t) which satisfies the Liouville-von 
Neumann equation 

ihp T (t) = [H(t),p T (t)}. (6) 

Both of the electronic leads and the thermal bath are 
assumed to be in equilibrium state all the time and de- 
scried by the time independent equilibrium density ma- 
trix pl and pb, respectively. Assuming the initial state 
as pr(0) = p(0)plPb, we can write the state at time t in 
the form pr(t) — p(t)pLPB under the Born approxima- 
tion. p(t) is the density matrix of the system composed 
of the two QDs and the mechanical oscillator. Iterating 
equation (6) in the interaction picture to the second or- 
der and performing trace over the leads {Ttl) and the 
bath (Ttb) variables, we obtain the master equation for 
the reduced density matrix of the system in the Markov 
approximation 

P" = Qo + Ql + Ql + Q12 + Ql (7) 

On the right hand side of equation (7), §q denotes evolu- 
tion term of the system in which QD2 is coupled to the 
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harmonic oscillator. q\ describes the contribution from 
direct tunneling by QD1 in the absence of QD2. gJj is 
just the right hand side of the master equation in our 
previous work Q (j . representing the contribution from 



vibration assisted transfer through QD2 alone. g\ 2 is co- 
herent term of the transport involving the two dots and 
g v d accounts for dissipation of the vibrational mode. Ex- 
plicit expressions of these terms read 



Qq = r^[eic|ci + fajjQa)a + (e 2 - hil(a^ + a))cjc 2 , p v ], (8) 
Ql = \ E [SS„(0, 0, mc\p< Cl - c lC \p v - p v Cl c\) + ES5(0, 0, l)(2c lP ^c t 1 - c\c x p* - p"c\ Cl )], (9) 

y=Lr 



1 c? ST V" (S» mi+ " 1+m2+ ™ 2 



Q2 = o e 



2 ^— ' ^— ' mi!ni!m2!n2! 

y— Z,r mi ,ni ,m2,ri2— 



+^22y( m i; n ii 2)(A m2n2 p"« + ^4„ imi p 1 ' ! ' A m2Tl2 — A^ n27l2 A nimi p v — p v A mini A m2Tl2 )], (10) 

& = V/ 2 y e «,*/a y iVL_ 

12 2 ^ ^ mini 

y—l,r m.n— 

x [2^(0, 0, l)(A+ nP ^ + Cl - p v c l A+ n ) + E$(0, 0, l)( C i/«A+„ - A+„c lP ") 

+S™ s (m I n, 2)(A+ m ^ + Cl - ^ nCl /) + S^(m, n,2)( Cl /a+ n - ^cMiJ] + ft.c, (11) 



and 

g v d = E m D[a^}p v + E out D[a]p v . (12) 

The degree of freedoms in the electronic leads and the 
thermal bath are assumed to be continuous with densi- 
ties of states Ny(^yk) and D(uj k ), respectively. The coef- 
ficients in the above equations corresponding to particle 
hopping into or out of the system are composed of inte- 
grals over these reservoir variables via 

E^Cmi.ni, = J d£, yk T l3y (£ yk )fy(t yk ) 

x<5(£ y fc - e* - (mi - ni)fuv ), 

x$(£ v k - e* - (mi - ni)fuJo), 
S m = J du} k ^(uj k )n B (uJ k )S(uj k - w ), 

and 

S ou * = J du k j(uj k )(l + n B (uJ k ))5(uj k -u) ). 



Here, T^y^yk) = 2irN y (£ yk )T* y T 3y , j(u k ) = 2nD(Lo k )g 2 
and i, j = 1,2. We have the Fermi-Dirac distribution 
function in lead y, f y (i yk ) = [e^y"~^/ kBT + l]^ 1 and 
the Bose-Einstein distribution function of the thermal 
bath, n B (u k ) = (e nUk/kBT - l)~ , where T is tem- 
perature and k B is the Boltzmann constant. In the 
above equation we have defined A~ m = C2(a') m (a) n , 
A tan = 4(at) m (a) n , where A+ n (A~ n ) describes that 
an electron hops into (out of) QD2 accompanied by cre- 
ation of m phonons and annihilation of n phonons. v, 
v+ = v + (1 + S y )/2 and v~ = v - (1 + S v )/2 indicate 
the number of electrons accumulated in the right lead. 
They are achieved by the following way: Trj j (dl k d rk pi J ) 

and TrL(d r kPLd rk ) contain different information about 
the number of electrons in the right lead when the num- 
ber is not infinite. Assuming v electrons are in the 
right lead, then the number of electrons in this lead 
can be expressed by p v f y (£, r k) = pTr L (dl k d rk p L ) and 
P v+1 fy{Crk) = pTr L (d rk p L d) rk ). In the same way, we 
have p v ~ 1 f y {^rk) — pTr L (dl k p L d rk ). The density matrix 
satisfies J2 v =o P" = P- The above method is equivalent 
to the counting approach in the many-body Schrodinger 
equation (33|, representing how many particles arrive at 
the collector. D[a]p v is a super operator acting on the 
density matrix p° by D[a]p v = ap v a* — (a)ap v +p v cva)/2. 
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3. Numerical treatment and current formula 

In this section, we give a brief introduction to our 
mathematical approach. We pay our attention to the 
character of stationary transport. Therefore, the elec- 
tron transmission can be conveniently described by the 
rate of electrons collected in the right lead. The system 
current is calculated according to the formula [34| 

oo 

I=-eJ2 v P V - ( 13 ) 

Here, P v = Tr mech [Tr char [p v }] is the probability of v 
electrons arrived at the right lead. The trace Tr mec h 
is taken over variables of the mechanical oscillator and 
Tr c har is taken over the degree of freedom of elec- 
tron occupation in the QDs. For the following numer- 
ical treatment, we consider the wide band approxima- 
tion and apply energy independent transmission rates 
Tijy — 2nN y T* y Tjy = 1,2 and y = l,r) and the 
damping rate 7 = 2irDg 2 . We assume r 12y and T 2 i y are 
real and satisfy Ti2 V — ^2iy = yTujT^j. The chemical 
potentials for the left and right electrodes are set to be 
/if = eV/2 and p r — —eV/2, respectively. 

The Fock state will be applied for the represen- 
tation of the master equation. And so the Hilbert 



Here, I\ is the current through QD1 alone, I2 is the cur- 
rent across the electromechanical junction in the absence 
of the reference arm. In fact, it is the same as the cur- 
rent directly derived from the master equation of the 
single- molecular junction [To| . And I12 is the interfer- 
ence part in terms of the off-diagonal density matrix for 
the electronic states. The values of density matrix el- 
ements are achieved by solving equation (7) under the 
condition p — 0. We project the equation in the basis of 
the Hilbert space as 

((m\(ji\(g + Qi + Q2 + Q12 + Qd)\ij)\n)) = 0, (18) 



space of the system is generated by the composite ba- 
sis {|00), |01), |10), |11)} ® {|0», |1», ..»), ...}, where ® 
means direct product. The state represents i elec- 
trons in QD1 and j electrons in QD2. \n)) is eigenstate 
of the nth excited level of the mechanical oscillator. In 
the Hilbert space, the system density matrix element is 
written as pijki, mn = ((to <g> (ji\p\kl) <g> \n))(i,j,k,l = 
0, l;m, n = 0,1,2,...). For any two given vibrational 
states to)), |n)) we have 16 density matrix elements 
PijM,mn(h j,k,l = 0,1) in terms of the electron states. 
However, just 6 of them are enough to describe the trans- 
port process since they can constitute a closed equation 
set for the system dynamics. These matrix elements are 

Pijij,mn and Pjiij : mn i^ij — 0, 

For the case of strong inter-dot Coulomb interaction, 
we assume that the state of two-electron occupation is 
not inside the transport window. In other words, the bias 
voltage is so low that only one electron passes the system 
at any time. As a consequence, the process involving 
the state pnn.mn is not contained in our equations [33j . 
Substituting equation (7) into equation (13), we reach 
the following expression for the current 

I = h+h+h2, (14) 

where 



(15) 



(16) 



(17) 

I 

((m\(ij\(go + gi + Q2 + Q12 + Qd)\ij)\n)) = 0, (19) 

where i, j — 0, 1 excluding the case i = j = 1. Then, one 
obtains a set of 5(N+ 1) 2 linear equations. N = 0, 1, 2, ... 
is the number of excited vibrational levels considered 
here. These equations can be solved by associating with 
the normalization condition ^2 V=Q P" = 1- For the nu- 
merical treatment, we take N = 18. It is a good ap- 
proximation for the low bias voltage, the weak dot-lead 
couplings and the finite oscillator damping rate applied 
here, as the contribution from the higher levels (n > N) 
of the vibrational mode is very small. 



h = eTr mech [E5#(0, 0, l)p 10 io - Eft r (0,0,l)poooo], 

I 2 = _ e Q V ^ 

2 ^— ' TOi!ni!TO2!ri2! 

mi ,ni ,m2 ,^2— 

xrr mech [S^(TO 1 ,n 1 ,2)(( a tH( a r i (at)^(a)™ 2 + (a*)™ 2 (a)™ V)" 1 (a)™ 1 V0101 
-S™ r (m 1 ,n 1 ,2)((at)«^(a) ro H« t r 2 (< 2 + (« t ) ro H«)" 2 (« t rH< 1 )/'oooo] ! 

Ii2 = y° 12 V ^^Tr mec/ Je^/ 2 (£^(0,0,1)(^ 
i — 4 mini 
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4. Results and discussions 



A. Phase relaxation and visibility 

In interference of two waves, visibility can be reduced 
not only by phase destruction of the waves but also due to 
difference of the absolute values of their amplitudes. The 
electromechanical systems substantially enhance electron 
transport for certain bias voltage [IMS]- Therefore, yo 
see the net contribution of phase relaxation to the in- 
terference fringe, we balance the amplitudes of waves 
in the two paths by taking the bare transmission rates 
of QD2 smaller than that of QD1. To this end, the 
bare tunneling rates for the reference path are set to be 
Tin — Tiir = O.Olw and for the target path taken as 
T121 — T*22r — 0.001312tJo- Then, the current in equa- 
tion (15) nearly equals to that in equation (16) with the 
small difference I\ — 1% < 10 _5 /ewo- In this case, we 
suppose that the absolute values of the two amplitudes 
corresponding to the two paths are almost the same. The 
current versus magnetic flux is plotted by the solid line 
in figure[2l It shows AB oscillation with the period of 2tt. 
Obviously, the interference current is not destructed to 
be vanished at its weakest points ((2n+l)7r, n is integer). 
It reveals that coherence of electron wave is influenced by 
the electromechanical vibration. Using the same way of 
balancing current amplitudes in the two paths, we give 
other three examples in figure [2] for different parameters. 
The low bias voltage eV = 3Hujq (red dotted line), high 
damping rate 7 — O.Suio (green dashed line) and small 
tunneling length a — x /X = 0.3 (blue dot-dashed line) 
weaken the effect from vibrational mode. As a result, the 
interference fluctuation is enhanced. 

Figure [2] shows the minimum and maximum values of 
the interference pattern is not shifted remarkably under 
different parameter variance. Using this feature, current 
visibility can be calculated easily by making a replace- 
ment as I max ~ I(4> = 0), I min ~ I(4> = 7r). It works 
under the conditions of £1 = £2 = and CI <C ujo- The 
visibility of interference fringe is given by the formula 
Visibility = (I ma x - Imin) / (Imax + Imin)- In figured 
(a), a substantial influence of the bias voltage to the in- 
terference visibility can be seen. For very low voltages 
eV < 2Hujq, there is no exited level contained in the 
transport window (eV), and visibility is close to unity. 
When applied voltage is close to zero, the correspond- 
ing current approaches to be vanished. It causes a small 
drop of visibility near the zero voltage. Excited states of 
the mechanical oscillator play an important role for the 
phase relaxation of electrons. Increasing the bias volt- 
age, excited levels of the vibrational mode are involved 
in the transport, which suppresses the visibility. At low 
voltage few discrete states of the vibration con- 

tribute to the transport and the visibility displays a step 
profile. The mechanical oscillation is naturally coupled to 
the thermal bath and it has an intrinsic life time which is 
the inverse of the damping rate 7. By rising the damping 
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FIG. 2: (Color on line) Current as a function of the magnetic 
flux through the AB ring. For the solid black line we take 
the transmission rates T221 ~ T22r = 0.001312ojo. For the red 
dotted line corresponding rates are T221 = ^22r ~ 0.002731o;o- 
For the green dashed line they are I^z = T22r = 0.003674o;o 
and for the blue dot-dashed line the transmission rates are 
T221 = T22r = 0.008091F The rest parameters are the same 
for all of the curves as ei = £2 = 0, Tm — Tn r — 0.01o;o, 
k B T = OmfkJo, xo/d = 0.003. 
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FIG. 3: (Color on line) (a) The visibility versus the bias volt- 
age (ei =£2=0, Tn; = Tn r = F 2 2l = r 2 2r = O.OlWo, 

a = 0.7, 7 = O.Olwo, k B T = 0.03Sw and x /d = 0.003). 
(b) The visibility as a function of the oscillator damping rate 

(eV = WlUJo, £i = £ 2 = 0, T 1U = Tllr = r 2 2i = T 2 2r = 

O.Olwo, a = 0.7, k B T = 0.03&j and x /d = 0.003). (c) The 
visibility as a function of the tunneling length (eV = 9&jo, 

£l = £2 = 0, Tin = Tn,. = T 2 2i = T 2 2r = O.OIcjo, 7 = O.Olwo, 

k B T = Omtkuo and x /d = 0.003). 



rate, the visibility increase can be observed as shown in 
figure [3] (b) . It is not hard to be understood that contri- 
bution from the mechanical motion would be decreased 
in the case of a high damping rate. The visibility is no 
longer enhanced obviously for damping rate 7 > 0.2u;o . 
It is in accord with the transition of the electromechani- 
cal system from so called shuttling regime into tunneling 
regime 0, 0]. The visibility is still not very high even 
at the quality factor ujq/j < 1. It implies, for the large 
damping rate, that coherence of electron does not obvi- 
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ously depend on the intrinsic lift time of the mechanical 
oscillator. In fact, the interference pattern is essentially 
affected by the strength of electron-phonon interaction 
which is determined by the parameter a — xq/X. In fig- 
ure [3] (c) , the visibility versus the coupling strength a 
is plotted. For a given oscillator with mass m and fre- 
quency ujq, the zero point uncertainty xq is fixed, and 
the coupling strength is mainly related to the tunneling 
length A. For infinite large tunneling length A -> oo we 
have a — > 0. In this case, Qi in equation (7) is close to 
the form of Q\ and the effect of vibration in g 2 and gi2 
is approached to disappear. As a consequence, tunneling 
between the two electrodes and QD2 is almost indepen- 
dent of the dot displacement. Therefore, we obtain the 
visibility close to be unity as illustrated in figure [3] (c) . 

In general, a large current induced by the vibrational 
junction is a reason of the reduced visibility in the AB 
ring. For instance, when one takes the same bare tunnel- 
ing rates for the two paths as shown in figure[31 the proba- 
bility of an electron passing the target arm is much larger 
than that of the electron propagating through the refer- 
ence arm. There is another probable reason for the weak 
interference, namely phase shift of electron waves, and 
this is also our main interest in the present paper. The 
propagating of an electron wave through the vibrational 
junction gives rise to many scattering excited states. Es- 
pecially at a sufficiently high applied voltage, the electron 
is in superposition of a large number of single-particle ex- 
cited modes associated with the electron-phonon interac- 
tion. These excited states are characterized by phase ac- 
quirement related to absorbtion and emission of phonons. 
One can expect that scattering in the space of positive 
phase shifts symmetrically happen in the space of nega- 
tive phase shifts. As a consequence, interference of all the 
scattering waves does not exhibits global phase shift be- 
tween the two paths. This property is valid for the same 
dot levels, E\ — £2 and the weak charge-field coupling Ml. 
In the next subsection let us discuss the case E\ 7^ £2. 
The influence of the charge-field coupling strength to the 
electron coherence has been previously considered in a 
similar system [3l], HI] . 



B. Coherent phase shift 

As we mentioned above, there is no global phase shift 
when an election propagating through the single-dot elec- 
tromechanical system. But it does not means there is no 
phase shift when one component of the electron wave 
transports through any individual level of the system. In 
order to observe the phase change of propagating electron 
wave through the target system, we change gate voltage 
in the reference arm and sec the variance of AB interfer- 
ence oscillation. As illustrated in figure 2J the pattern 
of the interference oscillation is shifting continuously in 
one direction when the resonant level of the reference 
arm is moving. The phase shift breaks original symme- 
try of the interference fringe for the replacement of </> by 




FIG. 4: (Color on line) AB interference oscillation as a func- 
tion of the resonant level of QD1 (eV — llftwo, £2 = 0, 
Tiii = Tn,. = T 2 2i = r 2 2r = O.Olwo, a = 0.7, 7 = O.Olwo, 
k B T = 0.03?k<;o and x /d = 0.003) 



— 4>. It is induced by detuning between the two QDs. 
In the AB interferometer of electron transport, interfer- 
ence is remarkably strong only when the energy level in 
one path is close to that in another path [35]. In other 
words, propagating waves in the two path are required 
to be oscillating in (at least nearly) the same frequency. 
Although the QD1 is detuned from the electronic level of 
the molecular junction, the molecular system still has en- 
ergy levels provided by the mechanical oscillator. There- 
fore, interference is not disappeared in the case of the 
detuning except some phase shift. 

The propagating wave in the reference arm only in- 
terferes with the wave in the molecular junction whose 
resonant energy (£2 + Ae) is the same as the resonant 
energy (£1) of the reference arm. Here, Ae is defined 
as the energy acquired or lost by an electron due to its 
inelastic scattering on the vibrating QD. Therefore, in 
figure 21 the phase shift corresponding to the detuning 
£1 — £2 represents phase change of the sub transmission 
amplitude whose energy is £1 = £2 + Ae in the molec- 
ular junction. The total amplitude of electronic wave 
transferring through the molecular junction is, of course, 
superposition of all the sub transmission amplitudes with 
different resonant energies. 

By choosing particular detunings E\ —£2 in figure[5j we 
analyze quantitative phase shift, especially for the reso- 
nant levels. Without loss of generality, the zero point 
energy is taken at £2 = 0. The numerical results in 
figure [5] show the phase shift AO of the transmission 
amplitude with energy £2 + A£ in the molecular junc- 
tion roughly satisfies the relation (£1 — S2)/Mjq — A9/ir. 
When Ae = AnhujQ (An = 0,1,2,...), the relation be- 
comes An ~ A9/tt. An is defined as the net number 
of phonons involved in the inelastic transport. It inti- 
mates that the phase difference of propagating waves cor- 
responding to two adjacent vibrational levels is w. This 
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off phase character is analogous to the phenomenon de- 
scribed by the Friedel sum rule [36| . The sum rule relates 
the phase shift of a scattering electron to the number of 
states in the energy interval due to the scattering. How- 
ever, the electron number accumulated in the impurity 
which is described by the general Friedel sum rule is re- 
placed by the phonon number involved in the electron 
transfer in our present system. 

The phase shift is very sharp when the energy of in- 
cident electron sweeping over the resonant levels [37| . 
However, in our model the tunneling is dependent on the 
position of mechanical oscillator, so the phase change is 
continuous and very smooth. It is investigated by the 
fact that the position dependent tunneling causes inelas- 
tic process, which improves decay of the oscillating QD 
and broadens the energy levels of the system [38| . 

From figure we know that the changes in the applied 
voltage, the electron-phonon coupling a and the damping 
rate of the oscillator do not induce global phase shift in 
the AB interferometer. Therefore, the definitive phase re- 
lation of tt difference between two adjacent levels is inde- 
pendent of these parameters so long as they are properly 
taken that the discrete levels of the mechanical vibration 
effectively contributes to the electron transport. For in- 
stance, on one hand the applied voltage should be large 
enough that at least one excited level of the oscillator is 
included in the transport window. On the other hand, 
the voltage is not too large so that the feature of discrete 
levels involved in the tunneling is obviously manifested. 
In fact, the phase shift is just related to the unit quanta 
of the mechanical oscillator as shown in figure O 

According to the above analysis, the neighboring res- 
onant levels in the molecular vibrational junction are off 
phase by tt. It is the character of one dimensional quan- 
tum system which is considered in our model. Since, 
in one dimensional system, the upper energy level has 1 
more wave function node than the lower one, and each 
node changes the phase of transmission amplitude by tt. 
This property may be not true if the system is not strictly 
one dimensional [39(. In the experiment of AB interfer- 
ometer where a fixed QD is embedded in one of the arms, 
the phase behaviors are the same for all resonant levels 
of the QD (l^4l6| . Namely, all the resonant levels are 
in-phase. It is different from present effect found in the 
electromechanical system, where all the vibrational levels 
are coherently correlated with definitive phase difference 
of tt. The phase shift varies from to any large value, 
depending on the net number of phonons involved in an 
electron tunneling. 

The reason of the visibility depression mentioned in 
the last subsection becomes more clear now. Actually, an 
electron takes all the channels of the discrete vibrational 
levels which are involved in the transport process. There- 
fore, interference not only occurs between the propagat- 
ing waves in the two paths, but also occurs among the 
waves taking different channels of the vibrational junc- 
tion. As we analyzed above, any two neighboring chan- 
nels have a phase difference of tt. The wave functions 




FIG. 5: Currents versus the variation of magnetic flux are 
plotted by changing the gate voltage of QD1 (eV = 9ftuo, 

£2 = 0, Tni = Tllr = T 2 2! = r 2 2r = O.OlCJO, « = 0.7, 

7 = 0.01w , k B T = Omhujo and x /d = 0.003) 



taking different vibrational levels destructively interfere 
because of the phase differences. It is the reason of phase 
relaxation in the AB interference due to the vibrational 
junction (see figure [2]). It has been shown in a double- 
QD two-electron AB interference that two components of 
conductance oscillations with the same amplitude cancel 
each other due to their phase difference of tt • Since 
two components of the conductance are the same in am- 
plitude, the final conductance disappears in their system. 
In the present case, electron occupation probabilities on 
the resonant levels of the electromechanical system are 
not the same. Therefore, there is net current remained 
in the system, but it is not fully coherent. In fact, the in- 
terference between the different channels is also reflected 
in the direct transmission of charge through the elec- 
tromechanical system [lol ]. The current calculated from 
the scheme considering both diagonal and off-diagonal 
couplings between the vibrational mode and the electron 
tunneling is remarkably lower than that obtained by the 
approach in which only diagonal terms are taken into ac- 
count. This current suppression is related to the destruc- 
tive interference between different transport channels. 

5. Conclusions 

Electrons propagating through the single-molecular vi- 
brational junction are dephased. It is caused by the elec- 
tron scattering on the excited levels of the vibrational 



8 



mode. However, the interference fringe of the AB in- 
terferometer is not absolutely destroyed by the nano- 
electromechanical system. The visibility is sensitive to 
the applied voltage, the oscillator damping rate and the 
tunneling length. The transmission amplitudes corre- 
sponding to channels of the vibrational resonant levels 
are coherently correlated via any neighboring channels 
have a definitive phase difference of it. Because of the 
phase shifts between the resonant levels in the electrome- 
chanical junction, different branches of the transmission 
waves destructively interfere with each other. As a con- 
sequence, the electron tunneling through the system ap- 
pears to be not fully coherent. The character of the phase 
difference of n is robust with respect to the wide range of 



the bias voltage, the tunneling length and the life time of 
the vibrational mode. It just depends on the frequency 
of the mechanical oscillator. This work would provide 
a guidance for the experimental observation of dephas- 
ing in electron transport through a vibrational molecular 
junction and phonon assisted conductors. 
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